PROCEEDINGS 

0F SCIENCE 



Thermodynamic properties of QCD in external 
magnetic fields 



CO 



G. S. Bali 1 , F. Bruckmann* 1 , M. Constantinou 2 , M. Costa 2 , G. Endrodi 11 , Z. Fodor 3 , 



S. D. Katz 4 , S. Krieg 3 5 , H. Panagopoulos 2 , A. Schafer 1 , K. K. Szabo 3 

We consider the effect of strong external electromagnetic fields on thermodynamic observables 
in QCD, through lattice simulations with 1 + 1 + 1 flavors of staggered quarks at physical quark 
^vq masses. Continuum extrapolated results are presented for the light quark condensates and for 

their tensor polarizations, as functions of the temperature and the magnetic field. We find the 
light condensates to undergo inverse magnetic catalysis in the transition region, in a manner that 
the transition temperature decreases with growing magnetic field. We also compare the results 
to other approaches and lattice simulations. Furthermore, we relate the tensor polarization to the 
spin part of the magnetic susceptibility of the QCD vacuum, and show that this contribution is 
diamagnetic. 
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1. Introduction 

Strong electromagnetic fields represent an important ingredient for the description of various 
strong interaction systems including non-central heavy-ion collisions, the interior of dense neutron 
stars (magnetars) and the cosmology of the early universe. The external magnetic field in these 
situations is expected to have a magnitude reaching up to the hadronic scale, such that its cou- 
pling to charged quarks and the color interaction between the latter compete with each other. This 
competition gives rise to an interplay between the strong and electromagnetic dynamics, producing 
various new phenomena. We will restrict ourselves to the idealized situation of constant external 
magnetic fields in thermal equilibrium. In fact, such magnetic fields probe the QCD vacuum in 
several aspects, by affecting its fundamental properties like chiral symmetry breaking and restora- 
tion, (de)confinement and hence the phase diagram, as well as the vacuum polarization and other 
quantities. 

At vanishing quark density, the transition of QCD at the physical values of the quark masses 
is a crossover [1] with transition temperatures, that may depend on the observable used for its def- 
inition. Here we focus on quark condensates, pseudo-order parameters related to the breaking of 
chiral symmetry, and the magnetic susceptiblities, relevant for various phenomenological applica- 
tions, in external magnetic fields. The renormalization of both will be defined, numerical results 
are presented for zero and nonzero temperatures and compared to results from other approaches. 
While the most important effects in the condensates are the inverse magnetic catalysis around the 
transition and the transition temperature decreasing with the magnetic field, the magnetic suscep- 
tibility contains information about the spin contribution to the para- or diamagnetic nature of the 
QCD vacuum. In these talks we summarize our publications [2, 3, 4], adding another plot for the 
continuum limit of the condensate (Fig. 4, right panel) and a few new comments. 

1.1 Magnetic catalysis and inverse magnetic catalysis 

A particularly pronounced effect of the magnetic field on the QCD dynamics is the so-called 
magnetic catalysis (MC) mechanism - the fact that the condensate increases with the external field 
B [5, 6] at zero temperature. This behavior has been observed in various effective theory and model 
calculations (for a recent review, see Ref. [7]), and has also been confirmed by lattice simulations 
in quenched theories [8, 9], at larger than physical pion masses in Nf = 2 QCD [10, 1 1] and in the 
N f = 4 SU(2) theory [12], and at physical quark masses in full N f = 2 + 1 QCD [2]. 

While all effective descriptions and lattice simulations agree about MC at T = 0, the situation 
becomes more complicated at finite temperature. As lattice simulations at physical quark masses 
with results extrapolated to the continuum limit show, the condensate exhibits a non-monotonic 
dependence on the magnetic field around T c , with a certain region, where it actually decreases 
with growing B [2]. This peculiar behavior we refer to as inverse magnetic catalysis (IMC). A 
consequence of this non-trivial dependence \j/y(B,T) is that the transition temperature T C (B) is 
reduced as the magnetic field increases. As it turns out, employing physical quark masses in the 
simulation and extrapolating the results to the continuum limit - as was done in our studies - is 
essential, as we will discuss. It is also highly important to address the differences between the lattice 
QCD results and model and chiral perturbation theory (#PT) predictions, especially since the latter 
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methods can be used to investigate regions that are not easily accessible to lattice simulations, e.g., 
QCD at a non-vanishing baryon density. 

1.2 Paramagnetism and diamagnetism 

Besides the influence of the magnetic field on chiral symmetry breaking, a fundamental aspect 
of the coupling between the strong dynamics and B is given in terms of the response of the free 
energy density, 

f = -^log&, (1.1) 



where 3f is the partition function of the system and V the (three-dimensional) volume. Due to 
rotational invariance, the B-dependence of / is to leading order quadratic, characterized by the 
magnetic susceptibility of the QCD vacuum, 



d 2 f 



d(eB) 2 



(1.2) 



eB=0 



which is a dimensionless quantity (here e > denotes the elementary charge). A positive suscep- 
tibility indicates a decrease in / due to the magnetic field, that is to say, a paramagnetic response. 
On the other hand, E, < is referred to as diamagnetism [13]. 

In the functional integral formalism of QCD, the susceptibility is readily split into spin- and 
orbital angular momentum-related terms [4], according to 



/ 



(1.3) 



with contributions from each quark flavor / with electric charge qf and mass mj. For a constant 
magnetic field B = Fxy in the positive z direction, 



g//g d 
2m f d(eB) 



->/JV 



eB=0 



(1.4) 



is given by an analogous expression with oVy replaced by a generalized angular momentum also 
present for spinless particles [4]. Eq. (1.4) constitutes an important relation which, to our know- 
ledge, has not been recognized previously in this context. Its derivation from the quark determinant 
and the corresponding Dirac operator is given in Ref . [4] . 

Eqs. (1.3) and (1.4) also show how the magnetic response of the free energy density is re- 
lated to the breaking of Lorentz symmetry by the external magnetic field. The presence of the 
preferred direction (B \\ z) induces a nonzero expectation value for the tensor polarization operator 
WfO^vW/, which appears in the spin contribution to the susceptibility, Eq. (1.4). To leading order, 
this expectation value is proportional to the field strength, and, thus, can be written as [14] 



(YfOxyWf) = q f B ■ (%\jrf) • Xf = Qf* ■ V, 



(1-5) 



where the expectation value of the quark condensate appeared. In the literature, %f i s referred to 
as the magnetic susceptibility of the condensate (for the quark flavor /). In what follows, we will 
also use the term "magnetic susceptibility". Again, we stress that it constitutes only one of the two 
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contributions to the total susceptibility. We also define the tensor coefficient T/ as the product of 
the condensate and the magnetic susceptibility. At finite quark masses, it is advantageous to work 
with Zf instead of %f for reasons related to renormalization (see Sec. 2). 

The magnetic susceptibility %f, in the context of QCD, was first introduced in Ref. [14]. It is 
a highly relevant quantity for experiments, since it appears in various effects [15, 16, 17, 18, 19]. 
In the past, the magnetic susceptibility has been calculated using different methods [20, 21, 22, 23, 
24, 25, 26, 27, 28, 29, 30, 31]. On the lattice, the numerical value of %f was determined recently 
in the quenched approximation of two- [32] and of three-color QCD [9], in both cases without 
renormalization. In this talk we report about new results regarding the susceptibility, at physical 
quark masses in full QCD [4]. 

An important remark regarding the magnetic susceptibility, which was not mentioned in Ref. [4] , 
is in order here. At vanishing temperature, electric charge renormalization ensures that the free en- 
ergy to €?(B 2 ) is given solely by the energy of the external field itself, B^/2, where B r is the 
renormalized magnetic field, see e.g. [33, 34]. Since this term is independent of the properties of 
the medium (i.e., of the QCD vacuum), we do not take it into account, resulting in a vanishing 
zero-temperature susceptibility, 

S(r = o) = o, (i.6) 

showing that at T = 0, the spin- and orbital momentum contributions are equal and of opposite 
sign. At T > this relation does not hold anymore, and the two terms become independent. Note 
that <^(fi 4 ) terms are already present in the free energy at T = 0, but these do not contribute to t,. 



2. Observables and renormalization 

To realize the external magnetic field on the lattice, we implement the continuum U(l) gauge 
field satisfying d x A y — d y A x = B, using space-dependent complex phases (for our implementa- 
tion, see [2]). This discretization satisfies periodic boundary conditions in the spatial directions, 
and ensures that the magnetic flux across the x — y plane is constant (it is also quantized). 

We consider three quark flavors u, d and s. Since the charges and masses of the quarks differ, 
we have to treat each flavor separately; q u = — 2q c i = —2q s . Furthermore, we assume m u = m c i^z m s . 
The partition function in the staggered formulation then reads, 

2f=(@Ue-f is s Y[ [detM(U,q f B,m f )] l/4 , (2.1) 

f=u,d,s 

with M{U ,qB,m) = 0(U,qB) + mt being the fermion matrix and /3 = 6/g 2 the gauge coupling. 
The exact form of the action we use and details of the simulation setup are given in Ref. [2] and in 
references there. 

In this formulation, the expectation value of the quark condensate for the flavor / can be 
written as 

, , T d log 3f T , i . . , 

^> = = W (TtM-^qfrmf)). (2.2) 

Likewise, the expectation value of the tensor Dirac structure of Eq. (1.5) reads, 

{WfO^Wf) = (TrM- l (U,q f B,m f )a^ v ) . (2.3) 
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In order to determine the continuum limit of the above observables, their renormalization has 
to be performed. The quark condensate (at finite mass) is subject to additive and multiplicative 
renormalization. All of these divergences - being independent of the temperature and of the mag- 
netic field - cancel [2] in the following combination [3], 

E M (S,T) = [WuAB,T) - ¥¥t<A0,0)] + 1, (2.4) 

where, to obtain a dimensionless quantity, we divided by the combination M 2 F 2 which contains 
the zero-field pion mass M % = 135 MeV and (the chiral limit of the) pion decay constant F = 
86 MeV [35]. This specific combination enters the Gell-Mann-Oakes-Renner relation, 

2m ud - ¥ Y(0,0)=MlF 2 + --- . (2.5) 

The factor +1 is included in Eq. (2.4), so that the chiral limit of the condensate is fixed to 1 at 
T = B = 0, and approaches as T — > °°. At nonzero quark mass, £ M ^ will still start from 1 at 
T = B = 0. 

We also define the change of the renormalized condensate due to the magnetic field as 

A£ M (B, T) = E M (B, T) - I M (0, T) = [WuA B > T ) ~ WuA°> r )] ■ ( 2 - 6 ) 

Note that the yh/A(0,0) term has canceled from this difference. Again, our normalization is such 
that the change of the condensate is measured in units of the chiral condensate at B = and T = 0. 
This normalization will be advantageous when comparing the lattice results to %PT and model 
predictions, which are usually given in units of v/i/^(0,0). 

The renormalization properties of the tensor polarization are somewhat different. As a calcu- 
lation in the free theory shows, an additive divergence of the form Qi-qfB- mflog(m 2 a 2 ) appears 
in {\j/fa^ v Yf) [4]. Moreover, there is a multiplicative divergence as well, which will be canceled 
by the tensor renormalization constant Zj. Both Qi and Zt are independent of T and B (and in 
mass-independent schemes of m/). In the free theory, the coefficient of the logarithmic divergence 
is calculated to be £ T (g = 0) = 3/(47r 2 ) [4]. 

Due to the divergence structure of the condensate and of the tensor polarization, it is advan- 
tageous to consider the tensor coefficient T/ defined in Eq. (1.5). We notice that the operator 
1 —mfd/dnif eliminates the logarithmic divergence, and, thus, can be used to define an observable 
with a finite continuum limit, 

/ n \ 

_div 

dm fj 



T /= Uf Zr^T/Zr-if . (2.7) 



At finite quark mass, this is one possible prescription to cancel the additive logarithmic term. It has 
the advantages that the chiral limit of Xf is left unaffected, and that, together with the logarithmic 
divergence, scheme-dependent finite terms also cancel in this difference [4], such that the scheme- 
and renormalization scale-dependence of Zf resides solely in Zj. The latter renormalization con- 
stant we calculate perturbatively, in the MS scheme at a renormalization scale of /I = 2 GeV. For 
the details of the perturbative calculation, see Ref. [4]. 
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Figure 1: QCD phase diagram in the magnetic field-temperature plane. Shown is the transition temperature 
T C (S) obtained from the renormalized light condensates, Eq. (2.4), averaged over the u and d quark. 



3. Results 



We start with our main result, the QCD phase diagram in the B — T plane, Fig. 1 . The transition 
temperature from the renormalized condensate of the light quarks decreases by up to 20 MeV for 
magnetic fields up to eB ~ 1 GeV 2 . A similar behavior is found in the chiral susceptibility and 
the strange quark number susceptibility as well [2]. We stress that these quantities are given in the 
continuum limit extrapolated from = 6, 8, 10 lattices using physical quark masses. 

The transition remains a crossover up to our largest magnetic fields as we verify by the finite 
volume scaling and the relative change in the T-dependence of the chiral susceptibility, Fig. 2. 




Figure 2: Strength of the transition analyzed through the light susceptibility % u on N t = 6 lattices, for details 
see [2]. Neither the volume scaling (left) nor the relative change in the peak shape (right) show evidence for 
a strengthening of the transition with the magnetic field, up to eB = 1 GeV 2 . 
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Figure 3: Magnetic catalysis at zero temperature. Shown is the change in the average of the light conden- 
sates, Eq. (2.6), due to the magnetic field using five lattice spacings and the extrapolation to the continuum 
limit (left) and then in comparison with the NJL model [36] and chiral perturbation theory [37, 38] in the 
low-fi region (right). 



We start our discussion of the light condensate at zero temperature. The two panels of Fig. 3 
show our lattice data and the extrapolation to the continuum limit, and a comparison to the NJL 
model and chiral perturbation theory. We confirm the increase of (YfWf) w i tri growing magnetic 
field, the well-known magnetic catalysis. The NJL and #PT predictions agree quantitatively with 
our data for eB < 0.3 GeV 2 and eB < 0.1 GeV 2 , respectively. The limitation of these approaches 
is not unexpected. The breakdown of the #PT prediction, for example, may be understood from 
the fact that the QCD dynamics is less and less dominated by pions as B grows, and charged p 
mesons play a non-negligible role, cf. [34, 39]. 




Figure 4: Inverse magnetic catalysis at the QCD transition. The light condensate as a function of the 
magnetic field for fixed temperatures (left) and as a function of T for fixed B (right) displays a non-monotonic 
behavior with MC at zero temperature turning into IMC around T c . The transition temperatures T C (B) of this 
observable are obtained from the inflection points in the right panel and hence can be seen to decrease with B. 
The right plot is the continuum extrapolation of Fig. 1 in [40] (where a slighty different normalization was 
used). 
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Figure 5: Influence of the quark mass on the dependence T C (B) of the (bare) light susceptibility on N t = 6 
lattices. For physical masses (left) the maximum moves to lower temperatures with the magnetic field, 
whereas in the system with three quarks at the strange quark mass (right) no change in position of the 
maximum is visible. This means that the effect of a decrease of T c with B is washed out by heavier-than- 
physical quark masses. 

For temperatures around the transition, we show the dependence of the light condensate on the 
magnetic field in Fig. 4. As advertised, the behavior is non-mono tonic. In the temperature range 
150-175 MeV the condensate decreases with the magnetic field, which is the inverse magnetic 
catalysis. It is especially this behavior that is missing in almost all effective theories of QCD, see 
e.g. the finite temperature comparison to a PNJL model and ^PT in Fig. 4 of [3]; for a recent 
attempt see [41]. From the right panel of Fig. 4, the decrease of the transition temperature with the 
magnetic field is visible. 

For the comparison to other lattice simulations of finite temperature gauge theories with mag- 
netic fields in SU (3) [10] and SU (2) [12], we first stress that our lattice results have been extrap- 
olated to the continuum limit. Another crucial difference to those simulations is the mass of the 
quarks, which is physical only in our simulations. The consequences of this we have investigated 
by artificially lifting the masses of the light quarks to the strange quark mass in our system, per- 
forming new simulations. As Fig. 5 shows by virtue of the light susceptibility, no decrease in T c 
with the magnetic field can be found anymore. The effect of IMC is washed out by the heavy 
quarks as well, cf. Fig. 5 lower panels in [2]. This mass dependence should, at least partially, be 
responsible for the different results obtained in other lattice simulations [10, 12]. 

With the mass-degenerate light quarks employed here, isospin symmetry is exact at vanishing 
field. At non-vanishing magnetic fields, however, this symmetry is broken by the quark charges 
—2e/3 vs. e/3. Therefore, we also consider the difference of the light condensates, A£„ — AZj = 
£„ — Erf. As Fig. 6 shows, this observable takes on nonzero expectation values and its temperature 
dependence is similar to that of the average light condensate. 
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Figure 6: Difference of the light condensates, AL U — A£j = E u — E f / signalling isospin breaking, as a function 
of temperatures for different magnetic fields, in analogy to Fig. 4 right panel. 



In the second part of this section, we present lattice data on the tensor polarization, again 
continuum extrapolated, and compute the magnetic susceptibility Xf an d the tensor coefficient Xf 
as the leading order proportionality factors in Eq. (1.5). 

First of all, Fig. 7 clearly reveals a linear growth of the tensor polarization (}PfOxyWf) w i tn the 
magnetic field. Already from this figure it is clear that both %f an d T/ are negative, and that the spin 
contribution to the QCD vacuum is, thus, of diamagnetic nature. This holds for all temperatures up 
to 200 MeV. In solid state physics, where one is interested in the nature of positive energy electrons, 
the spin is associated with paramagnetism. Here, however, we analyze the response of the QCD 
vacuum, which is different, because, e.g., the vacuum energy for fermions is negative. For reviews 
about the subject we refer the reader to Refs. [42, 43, 44]. 

The corresponding tensor coefficients at zero temperature are plotted in Fig. 8 as a function 
of the light quark masses, that we varied between the physical value and the strange quark mass, 
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Figure 7: Negative tensor polarization as a function of the magnetic field for three temperatures near the 
transition. The agreement between the data for different spatial volumes at T = 141 MeV indicates negligible 
finite volume effects. 
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Figure 8: Mass dependence of the combination — ZjX,, in the MS scheme at renormalization scale jU = 2 
GeV. The coefficient of the logarithmic divergence is determined by fitting the data by Eq. (3.1) (solid lines). 



again through new dynamical simulations. These data have already been multiplied by the tensor 
renormalization constant Zj (and —1), but still contain a divergence of the form log (injd) as we 
have discussed before Eq. (2.7). Correspondingly, we fit these data with coefficients c using 

c + c l R + c 1 R\og{R 2 a 2 ), R = m ud /mfj\ c, = cf ] + c|V , (3.1) 

such that 

T' f = Cfo + 2c f2 m f . (3.2) 

We remark that the coefficients of the logarithms, c^/nP^P = 0.055(5) and c^/nt^ s = 0.072(6) 
are quite close to the free-field value of 3/(4n 2 ). As a result of these fits, we get the following 
values for the renormalized tensor coefficients in the chiral limit, 

< = -40.3(1. 4) MeV, x r d = -38.9(1.5)MeV, (3.3) 

and at physical quark masses, 

x r u = —40.7(1.3) MeV, x r d = —39 .4(1.4) MeV, < = -53.0(7.2) MeV. (3.4) 

In order to extract the magnetic susceptibilities in the MS scheme at pL = 2 GeV, we use a 
recent lattice determination of = (269(2) MeV) 3 [45] and (\jr s yr s ) / (WiYi) = 0.8(3) from 

[46]. This yields 

Xu = - (2.08 ± 0.08) GeV" 2 , % d = - (2.02 ± 0.09) GeV" 2 , % s = - (3.4 ± 1 .4) GeV~ 2 . 

(3.5) 

These values are in good agreement with the QCD sum rule calculations summarized and 
updated in [22], Xud = —2.1 1 (23) GeV -2 . Quenched lattice simulations, on the other hand, gave the 
unrenormalized values % ud = -1. 547(6) GeV~ 2 for two colors [32] and Xud = -4.24(18) GeV -2 
for three colors [9]. 

Finally, we show the corresponding temperature dependence of % r u in Fig. 9. Since the diver- 
gent term T dlv of Eq. (2.7) is independent of T, we subtracted it in the same way as at T = 0. Then 
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Figure 9: Temperature dependence of the negative renormalized light tensor coefficient in the MS scheme 
at a renormalization scale jj, = 2 GeV for three lattice spacings and the continuum extrapolation. 

we performed a combined fit of the form N t ~ 2 for different lattice spacings. We found the light ten- 
sor coefficient to decrease in the transition region like a quasi-order parameter. The corresponding 
transition temperature, T c = 162(3) (3) MeV, is similar to that of the chiral condensate. 

4. Summary 

We have presented results on equilibrium thermodynamics of QCD in external magnetic fields, 
obtained through lattice simulations. We have found the well-known magnetic catalysis at zero 
temperature to turn into the opposite - inverse magnetic catalysis - at temperatures around the 
transition. We have argued that performing the continuum limit and using physical quark masses 
are essential for obtaining these results. The inverse magnetic catalysis behavior still awaits confir- 
mation from other first principle methods, for instance from lattice simulations with other fermion 
discretizations. This also applies to the associated decrease of the transition temperatures with the 
magnetic field, which could be of phenomenological relavance for, e.g., heavy-ion collisions. Both 
phenomena seem to represent a challenge for effective theories of QCD. 

As a second key aspect we have calculated one of the Lorentz symmetry breaking observables 
induced by the magnetic field, the tensor polarization. It gives a diamagnetic contribution to the 
spin-related part of the magnetic susceptibility of the QCD vacuum. 
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